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■ ABSTRACT 

Aims. Based on a new set of sulphur abundances in very metal-poor stars and an improved analysis of previous data, 
we aim at resolving current discrepancies on the trend of S/Fe vs. Fe/H and thereby gain better insight into the 

■ nucleosynthesis of sulphur. The trends of Zn/Fe and S/Zn will also be studied. 

Methods. High resolution VLT/UVES spectra of 40 main-sequence stars with —3.3 < [Fe/H] < —1.0 are used to derive 
S abundances from the weak A8694.6 Si line and the stronger AA9212. 9, 9237.5 pair of Si lines. For one star, the S 
abundance is also derived from the S I triplet at 1.046 /im recently observed with the VLT infrared echelle spectrograph 
^ ' CRIRES. Fe and Zn abundances are derived from lines in the blue part of the UVES spectra, and effective temperatures 

Q\ ' are obtained from the profile of the H/3 line. 

Results. Comparison of sulphur abundances from the weak and strong S I lines provides important constraints on non- 
\£) ■ LTE effects. The high sulphur abundances reported by others for some metal-poor stars are not confirmed; instead, 

£SJ ' when taking non-LTE corrections into account, the Galactic halo stars distribute around a plateau at [S/Fe]~ +0.2 dex 

with a scatter of 0.07 dex only. [Zn/Fe] is close to zero for metallicities in the range —2.0 < [Fe/H] < —1.0 but increases 
£ — \ to a level of [Zn/Fe]~ +0.1 to +0.2 dex in the range -2.7 < [Fe/H] < -2.0. At still lower metallicities [Zn/Fe] rises 

O ■ steeply to a value of [Zn/Fe]~ +0.5 dex at [Fe/H] = -3.2. 

Conclusions. The trend of S/Fe vs. Fe/H corresponds to the trends of Mg/Fe, Si/Fe, and Ca/Fe and indicates that 
' sulphur in Galactic halo stars has been made by a-capture processes in massive SNe. The observed scatter in S/Fe is 

SpH, much smaller than predicted from current stochastic models of the chemical evolution of the early Galaxy, suggesting 

that either the models or the calculated yields of massive SNe should be revised. We also examine the behaviour of 
5_l ■ S/Zn and find that departures from the solar ratio are significantly reduced at all metallicities if non-LTE corrections 

to the abundances of these two elements are adopted. This effect, if confirmed, would reduce the usefulness of the S/Zn 
ratio as a diagnostic of past star-formation activity, but would bring closer together the values measured in damped 
Lyman-alpha systems and in Galactic stars. 
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1. Introduction Other investigations point, however, to an increasing 

trend of [S/Fe] towards lower metallicities (Israelian & 
Despite several recent papers on sulphur abundances in Rebolo 2001; Takada-Hidai et al. 2002) with [S/Fe] report- 
Galactic stars, there is still no agreement on the trend edly reaching as high as +0.8 dex at [Fe/H] ~ -2. As dis- 
of [S/Fe] vs. [Fe/H] for halo stars. Some studies (Ryde & cussed by Israelian & Rebolo, such high values of [S/Fe] 
Lambert 2004, Nissen et al. 2004) indicate that [S/Fe] is may be explained if SNe with very large explosion energies 
approximately constant at a level of +0.3 dex in the metal- of E = (10 - 100) x 10 51 ergs (so-called hypernovae) make a 
licity range -3 < [Fe/H] < -1. Such a plateau-like over- substantial contribution to the nucleosynthesis of elements 
abundance of S with respect to Fe is also predicted from in the early Galaxy (Nakamura et al. 2001). An alternative 
Galactic chemical evolution models with Type II SNe as the explanation of an increasing trend of [S/Fe] towards lower 
dominant source of element production up to [Fe/H] - -1 metallicities has been proposed by Ramaty et al. (2000): 
(Chiappini et al. 1999, Goswami & Prantzos 2000); the de- assuming a short mixing time (~ 1 Myr) for supernovae- 
cline of [S/Fe] for galactic disk stars (Chen et al. 2002, Ryde synthesized volatile elements like oxygen and sulphur and 
2006) is then due to the release of iron from Type la SNe a longer mixing time (~ 30 Myr) for refractory elements 
at [Fe/H] > — 1. like Fe, high values of [O/Fe] and [S/Fe] are expected in 

the early Galaxy. 



* Based on observations collected at the European Southern In a recent P a P er h Y Caffau et aL ( 2005 )> more Puzzling 

Observatory at Paranal, Chile (programmes No. 67.D-0106, data on [S/Fe] m halo stars have been obtained. Both high 

73.D-0024 and CRIRES science verification program 60.A-9072) [S/Fe] ~ +0.8 dex and low [S/Fe] ~ +0.3 dex are found in 

** Table 1 and Appendices A, B, and C are only available in the metallicity range -2.2 < [Fe/H] < -1.0 (see Fig. 10 of 

electronic form at http://www.aanda.org Caffau et al.), suggesting a dichotomy of [S/Fe]. If real, this 
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points to a very complicated chemical evolution of sulphur 
in the early Galaxy. 

Additional problems have been revealed by Takeda et 
al. (2005), who find that S abundances derived from the 
weak A8694.6 S I line are systematically higher than S abun- 
dances derived from the stronger AA9212. 9, 9237.5 pair of 
lines when their new non-LTE corrections are applied. 

The uncertainty about the trend of [S/Fe] calls for fur- 
ther studies of sulphur abundances in halo stars. In this 
paper we present new sulphur abundances for 12 halo stars 
with —3.3 < [Fe/H] < —1.9 based on high-quality, near- 
IR spectra obtained with the VLT/UVES spectrograph. 
We have also re-analyzed UVES spectra of 28 stars with 
[Fe/H] < -1 from Nissen et al. (2004, hereafter Paper I) 
determining T e g- from the H/3 line in the same way as for 
the new stars and improving the derivation of S abundances 
from the AA9212.9, 9237.5 lines by taking into account the 
opacity contribution from the wings of the Paschen-zeta 
hydrogen line at 9229 A. In addition, possible non-LTE ef- 
fects are investigated by comparing S abundances obtained 
from the weak Si line at 8694.6 A with data from the 
AA9212.9, 9237.5 Si lines. 

We have also derived zinc abundances from the 
AA4722. 2, 4810.5 Zni lines. Both S and Zn are among the 
few elements which are not readily depleted onto dust in 
the interstellar medium of the Milky Way and are present 
in the gas-phase in near-solar proportions. For this rea- 
son, they are key to studies of metal enrichment in distant 
galaxies, particularly those detected as damped Lyman- 
alpha systems (DLAs) in the spectra of high redshift QSOs. 
Assuming that sulphur behaves like other ev-capture ele- 
ments and that Zn follows Fe, the S/Zn ratio may be used 
to date the star formation process in DLAs (Wolfe et al. 
2005). A clarification of the trends of both S and Zn for 
Galactic stars is important to test if these assumptions are 
correct. 

2. Observations and data reduction 

The 12 new programme stars were selected from the list of 
very metal-poor turnoff stars in Ryan et al. (1999) and 
observed with the VLT/UVES spectrograph (Dekker et 
al. 2000) in service mode during ESO period 73 (April - 
September, 2004). In addition to the programme stars, a 
fast-rotating bright early-type B star (either HR5488 or 
HR 6788) was observed on each night in order to be able to 
remove telluric H 2 lines in the 9212 - 9238 A region. 

The UVES setup was the same as described in Paper I 
except that no image slicer was applied when imaging the 
star onto the 0.7 arcsec wide entrance slit. Briefly, we men- 
tion that the dichroic mode of UVES was used to cover 
the spectral region 3750 - 5000 A in the blue arm and 
6700 - 10500 A in the red arm, in both cases with a spec- 
tral resolving power of A/ A A ~ 60 000. The blue region 
contains a number of Fen lines suitable for determining 
the iron abundance, and two Zni lines. In the red region 
we find the weak A8694.6 S I line as well as the stronger S I 
lines at 9212.9 and 9237.5 A. Typical S/N ratios are 300 
in the blue spectral region, 250-300 at the A8694.6 S I line, 
and 150-200 in the 9212 - 9238 A region. 

The spectra were reduced by using standard IRAF rou- 
tines for order definition, background subtraction, flat-field 
correction, order extraction with sky subtraction and wave- 
length calibration. Continuum fitting was performed with 




9210 9215 9220 9225 9230 9235 9240 
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Fig. 1. The VLT/UVES spectrum of HD 84937 in the spec- 
tral region 9210 - 9240 A overlaid with the B2III spectrum of 
HR 5488. Below is the spectrum of HD 84937 (shifted ~ 0.5 units 
in relative flux) after removal of the telluric lines using the IRAF 
task telluric to divide with a scaled version of the spectrum 
of HR5488. 



the IRAF task continuum using spline functions with a 
scale of a few A. Since many continuum windows are avail- 
able for these metal-poor stars, even in the blue part of the 
spectrum, this method works well. 

The quality of the 2004 spectra is about the same as 
shown in Fig. 2 of Paper I. As an additional illustration 
of the removal of telluric lines, we show the spectrum of 
HD 84937 in Fig. 1. This bright star was observed with 
particularly high S/N. As seen from the figure, the A9237.5 
S I line is completely overlapped by a telluric line, but after 
division with the scaled B-type spectrum of HR5488, the 
line emerges very clearly. Due to the radial velocity differ- 
ence of the two stars, the PC Hi line at 9229A is seen as 
a broad feature in 'emission' from HR 5488 overlaid with a 
narrower absorption line from HD 84937. The wings of the 
P^ line are lost and therefore the equivalent widths of the 
S I lines can only be measured relative to the local contin- 
uum. Hence, when using their equivalent widths to derive 
S abundances, the opacity contribution from the P£ line 
should be taken into account. 

Interference fringing is another problem in the near-IR 
part of the spectrum. The fringes on the UVES MIT CCD 
have an amplitude of 20 - 40% at 9000 A. After flat-fielding, 
residual fringes at a level of ±0.5% remain as seen from 
Fig. 2. As discussed by Korn & Ryde (2005), this makes 
it difficult to determine precise S abundances of halo stars 
from the weak A8694.6 Si line. B-type stars may, however, 
also be used to correct for the residual fringing as illustrated 
in Fig. 2. After division by the spectrum of HR5488, the 
spectrum of HD 84937 has a S/N ~ 500. Without such a 
correction, the error of the S abundance derived from the 
A8694.6 line can be large, especially for stars with [Fe/H] < 
-2. 

Equivalent widths of S, Fe and Zn lines were measured 
by Gaussian fitting or direct integration if the fit was poor 
and are given in Table A.l for the 12 stars observed in 2004. 
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Fig. 2. The VLT/UVES spectrum of HD 84937 in the spec- 
tral region 8685 - 8705 A overlaid with the B2III spectrum 
of HR5488. Below is the spectrum of HD 84937 (shifted down 
0.04 units in relative flux) after division with the spectrum of 
HR5488. 

The corresponding table for the Paper I stars is available at 
the CDS 1 . 
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Fig. 3. Observed H/3 profiles in the spectra of CD-24°17504 
(T cff = 6340 K, upper jagged line) and LP 815-43 (T cB = 6480 K, 
lower jagged line) compared to synthetic profiles (thin lines) 
computed for model atmospheres with logg = 4.3, [Fe/H] = 
-2.8 and T cff = 6000, 6300, 6600 K, respectively. An index, 
/3(UVES), measuring the strength of the H/3 line is defined from 
the total flux in the C bands relative to the total flux in the L 
bands (see Appendix B). 



3. Model atmospheres and stellar parameters 

The determination of the effective temperature {T e g) as 
well as the abundances of S, Fe and Zn is based on a- 
element enhanced ([a/Fc] = +0.4, a = O, Ne, Mg, Si, S, Ca, 
and Ti) ID model atmospheres computed with the Uppsala 
MARCS code. Updated continuous opacities (Asplund et 
al. 1997) including UV line blanketing are used. LTE is 
assumed both in constructing the models and in deriving 
T e g and the abundances. Convection is treated in the ap- 
proximation of Henyey et al. (1965) with a mixing- length 
parameter of omlt = l/H p = 1.5, and a temperature distri- 
bution in the convective elements determined by the diffu- 
sion equation (parameter y — 3/4 ir 2 — 0.076). The reader 
is referred to Ludwig et al. (1999) for a summary of the 
parameters used in different versions of the mixing length 
theory. 

3.1. Effective temperature 

Instead of deriving T c g from the (b — y) and (V — K) colour 
indices as in Paper I, we decided to use the profile of the H/3 
line in the blue part of the UVES spectra. The advantage 
is that the H/3-based T e g is not affected by possible errors 
in the interstellar reddening; in general the new stars are 
more distant than the stars in Paper I and hence potentially 
more affected by interstellar absorption. 

The H/3 line is well centered in an echelle order having 
a width of approximately 80 A. After flat-fielding, the con- 
tinua of the two adjacent echelle orders were fitted by low 
order spline functions and the continuum of the H/3 order 
was determined by interpolation between these functions 
in pixel space. The H/3 echelle order is not wide enough 

1 http : //cdsweb . u-strasbg . f r/cgi-bin/qcat?/A+A/415/993 



to reach the true continuum, but the method allows one 
to use the inner ±30 A of the profile to derive T e g from a 
comparison with synthetic H/3 profiles for a grid of model 
atmospheres. In addition to the 12 new stars, T e g was also 
re-determined in this way for stars from Paper I. 

To illustrate this method, Fig. 3 shows the H/3 pro- 
files in the spectra of two stars compared to synthetic pro- 
files calculated as described in Barklem et al. (2002) using 
Stark broadening from Stehle & Hutcheon (1999) and self- 
broadening from Barklem et al. (2000a). A few metallic 
absorption lines have been removed by interpolating the 
H/3 profile across these lines. An index, /3(UVES), defined 
as the ratio of the total flux in the two C-bands to the to- 
tal flux in the two L bands is used to determine T c g (see 
Appendix B). Hence, we are not using the position of the 
true continuum in determining T e g. The center of the H/3 
line is also avoided, because the fit between theoretical and 
observed profiles is poor around the line center due to non- 
LTE effects (Przybilla & Butler 2004). 

For stars with [Fe/H] > —1, the metal lines in the H/3 
line are so strong that the removal of these lines by inter- 
polation of the profile across the lines is unreliable. Hence, 
six stars with [Fe/H] > — 1 from Paper I are not included 
in the present paper. 

The symmetric appearance of the observed profiles 
shows that the procedure of rectifying the H/3 echelle order 
has worked well. In fact, the difference between T e g deter- 
mined from the left and the right part of the H/3 profile 
never exceeded 30 K. The internal stability of the method 
is also very good; the effective temperatures determined 
for a given star observed on different nights agreed within 
±20 K. Another advantage is that the profile of the H/3 line 
varies very little with gravity and metallicity. Hence, the 
errors in these parameters have only a small effect on the 
derived T e g. For a group of stars with similar atmospheric 
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Fig. 4. The difference of effective temperatures derived from 
(V — K)o and H/3 using three different calibrations of (V — K)o. 
The comparison includes only 38 stars; LP 65f-4 falls outside 
the range in (V — K)o where the calibrations are valid, and 
CD — 7f°f234 is missing Stromgren photometry; hence, its inter- 
stellar reddening excess, E(V — K), could not be derived. Three 
stars known to be single- lined spectroscopic binaries are marked. 



parameters, like metal-poor turnoff stars, one might there- 
fore expect that differential values of T e g can be determined 
with a precision of about ±30 K from the H/3 profile. 

The systematic error of T c g is, however, larger. The fit 
between the theoretical and observed profiles is not per- 
fect in the central region of the H/3 line, and the estimated 
value of T c ff therefore depends to some extent on where 
the L-bands of the H/3 index are placed. The temperature 
structure of the model atmospheres also plays a role. If the 
mixing- length parameter is decreased from l/H p = 1.5 to 
1.0, the derived T e g- changes by about — 25 K for metal- 
poor stars at the turnoff (Toff ~ 6400 K) and about — 60 K 
for cooler stars (T cff - 5800 K). 

The derived effective temperatures are given in Table 
1. Ten of the stars are in common with Asplund et al. 
(2006), who determined T c ff from the Ha profile in inde- 
pendent UVES spectra. The average difference, T e g (H/3) 
— T c ff (Ha) is 64 K with a rms scatter of the deviation of 
±28 K. This small scatter confirms that precise differen- 
tial values of T e ff can be determined from Ha and H/3 for 
stars with similar atmospheric parameters, but the system- 
atic difference indicates that the absolute values of T e g are 
more uncertain. Interestingly, the systematic difference of 
T e ff (H/3) and T c g (Ha) vanishes if the mixing-length pa- 
rameter is decreased to l/H p = 0.5. The same conclusion 
was reached by Fuhrmann et al. (1993) for a set of Kurucz 
model atmospheres with the convective structure param- 
eter y — 0.5. Given that Barklem et al. (2002) were not 
able to obtain a satisfactory fit to the Ha and H/3 line 
profiles in the solar spectrum for any combination of the 



convection parameters l/H p and y, it seems, however, pre- 
mature to change these parameters from the standard val- 
ues of the MARCS models. Clearly, the model atmospheres 
and/or the line broadening theory need to be improved. 
Although it is often argued in the literature that the wings 
of Balmcr lines are formed in LTE for late- type stars, very 
recent work suggests that this may in fact not necessarily 
be true (Barklem 2007). Unfortunately, the degree of de- 
partures from LTE depends critically on still poorly known 
inelastic collisions with other neutral H atoms. Some exist- 
ing atomic calculations imply LTE while others yield sig- 
nificant non-LTE effects. A full 3D non-LTE analysis of 
Balmer line profiles in cool stars with improved collisional 
data is urgently needed to resolve this potentially signifi- 
cant systematic error. 

The problem of systematic errors in the effective tem- 
peratures of metal-poor stars is also evident when com- 
paring T c ff (H/3) with effective temperatures derived from 
the (V — K) index using three different colour- T e g calibra- 
tions (see Fig. 4). Here, the V magnitude was taken from 
Stromgren photometry (see Table B.l), and K from the 
2MASS catalogue (Skrutski et al. 2006). Correction for red- 
dening was applied according to the relation E(V — K) = 
3.8 E(b-y) (Savage & Mathis 1979), where E(b-y) is de- 
rived as described in Appendix B. 

The (V — K)o calibrations of Ramirez & Melendez 
(2005b) and Masana et al. (2006) refer directly to the 
2MASS K s magnitude system, whereas the Alonso et al. 
(1996) calibration refer to the original Johnson (1966) Kj 
system. Hence, before applying the Alonso et al. calibra- 
tion, we converted K s to Kj using equations in Alonso et 
al. (1994) relating Kj to K C it (Elias et al. 1982) and the 
transformation from K s to Kcit derived by Cutri et al. 
(2006). These transformations all contain a small (J—K) 
colour term, but since our sample of stars are confined to 
a narrow range in J — K , we get in a good approximation 
(V — Kj) ~ (V — K s ) — 0.035 for metal-poor stars in the 
turnoff region. If this transformation had not been applied, 
the effective temperatures derived from (V — K) with the 
Alonso et al. calibration would have been about 65 K lower. 

As seen from Fig. 4, there is good agreement between 
T e ff from H/3 and from (V — K)q when using the Alonso 
et al. (1996) calibration. Excluding the three known binary 
stars, the mean difference, T e ff (V — K)o — T e g (H/3) is — 11± 
72 K. The scatter agrees well with the error of T c ff (V — K )q 
expected from errors in the reddening correction. Still, there 
is a tendency for the residuals to increase with decreas- 
ing [Fe/H]. This tendency is even more pronounced for the 
(V — K)o calibrations of Ramirez & Melendez (2005b) and 
Masana et al. (2006). When using the (V— K)q calibration 
of Ramirez & Melendez, stars with [Fe/H] > —2.0 have a 
mean deviation T e g (V—K)o — T e g (H/3) of about — 50 K, 
whereas stars with [Fe/H] < —2.5 have a mean deviation 
of +124K. Hence, our H/3-based effective temperatures do 
not confirm the hot T c ff scale of very metal-poor turnoff 
stars derived by Ramirez & Melendez (2005a) from their 
application of the infrared flux method. 

Four of the stars in Table 1 are known to be single- 
lined spectroscopic binaries, i.e. G 66-30 (Carney et al. 
2001), G 126-62 and G59-27 (Latham et al. 2002), and 
CD-71°1234 (Ryan et al. 1999). As discussed by Ramirez 
et al. (2006) in the case of G 126-62, these stars may have a 
cool companion that causes too low a value of T e ff to be de- 
rived from {V-K) . Indeed, G 126-62 and G 59-27 exhibit 
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Table 1. Derived stellar parameters and abundances of sulphur and zinc. LTE sulphur abundances are given for each of the three 
S I lines. The LTE Zn abundance is the average of the values derived from the AA4722.2, 4810.5 Zn I lines. In calculating [S/Fe] and 
[Zn/Fe], solar abundances loge(Fe)o = 7.50, loge(S)o = 7.20 and loge(Zn)0 = 4.60 have been adopted. The non-LTE values of 
[S/Fe] and [Zn/Fe] are based on the non-LTE corrections calculated by Takeda et al. (2005) for a hydrogen collisional parameter 
S H = 1. 
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1.5 
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0.16 
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-0.01 


0.05 
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1.5 


4.85 


4.85 




0.19 


0.07 


2.25 


0.19 


0.26 
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1.2 
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6157 


4.08 


-1.59 


1.5 


6.03 


5.92 


5.74 


0.29 


0.21 


3.02 


0.01 


0.07 


HD 140283 


5849 


3.72 


-2.38 


1.5 


5.07 


5.06 


5.17 


0.28 


0.20 


2.30 


0.08 


0.14 


HD 160617 


6047 


3.84 


-1.75 


1.5 


5.81 


5.80 


5.85 


0.37 


0.28 


2.82 


-0.03 


0.04 


HD 179626 


5881 


4.02 


-1.12 


1.4 


6.37 




6.31 b) 


0.26 


0.17 


3.48 


0.00 


0.04 


HD 181743 


6044 


4.39 


-1.87 


1.5 


5.62 


5.62 




0.29 


0.23 


2.81 


0.08 


0.12 


HD 188031 


6234 


4.16 


-1.72 


1.5 


5.76 


5.75 


5.78 


0.28 


0.21 


2.91 


0.03 


0.09 


HD 193901 


5699 


4.42 


-1.10 


1.2 


6.29 


6.30 


6.34 


0.21 


0.17 


3.38 


-0.12 


-0.10 


HD 194598 


6020 


4.30 


-1.15 


1.4 


6.26 


6.30 


6.20 b) 


0.20 


0.13 


3.37 


-0.08 


-0.04 


HD 215801 


6071 


3.83 


-2.28 


1.5 


5.28 


5.16 




0.30 


0.17 


2.41 


0.09 


0.17 


LP 815-43 


6483 


4.21 


-2.71 


1.5 


4.87 


4.81 




0.35 


0.17 


2.17 


0.28 


0.38 
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CD -24° 17504 


6338 


4.32 


-3.21 


1.5 


4.51 






0.52 


0.30 


1.96 


0.57 


0.66 


CD-71°1234 a) 


6325 


4.18 


-2.38 


1.5 


5.13 


5.19 




0.34 


0.22 


2.34 


0.12 


0.20 


CS 22943-095 


6349 


4.18 


-2.24 


1.5 


5.36 


5.33 




0.38 


0.26 


2.49 


0.13 


0.20 


G 04-37 


6308 


4.25 


-2.45 


1.5 


5.05 


4.98 




0.26 


0.15 


2.29 


0.14 


0.21 


G 48-29 


6482 


4.25 


-2.60 


1.5 


4.94 


5.03 




0.39 


0.23 


2.13 


0.13 


0.23 
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4.23 


-1.93 


1.5 


5.57 


5.57 




0.30 


0.20 


2.73 


0.06 


0.12 


G 126-52 


6396 


4.20 


-2.21 


1.5 


5.21 






0.22 


0.09 


2.50 


0.11 


0.19 


G 166-54 


6407 


4.28 


-2.58 


1.5 


5.04 


5.06 




0.43 


0.29 


2.10 


0.08 


0.16 


HD 84937 


6357 


4.07 


-2.11 


1.5 


5.47 


5.38 


5.40 


0.33 


0.23 


2.55 


0.06 


0.15 


HD 338529 


6373 


4.03 


-2.26 


1.5 


5.36 


5.32 




0.40 


0.25 


2.48 


0.14 


0.23 


LP 635-14 


6367 


4.11 


-2.39 


1.5 


5.12 


5.07 




0.29 


0.14 


2.36 


0.15 


0.24 


LP 651-4 


6371 


4.20 


-2.63 


1.5 


4.83 






0.26 


0.11 


2.04 


0.07 


0.16 



Single-lined spectroscopic binary star. 

The weaker A8694.0 S I line is included in the determination of this S abundance. 



some of the largest deviations in Fig. 4, whereas G 66-30 
does not show a significant deviation. This star has an un- 
usually high effective temperature (T c s— 6470 K) for a halo 
star with [Fe/H]= —1.48, and is classified as a blue straggler 
by Carney et al. (2001). 



3.2. Surface gravity 

As in Paper I, the surface gravity is determined from the 
fundamental relation 

log A = log J^L + 4 log + 0.4(M bo] - M bol;0 ) (1) 

where M is the mass of the star and Mboi the absolute 
bolometric magnitude. 

The Stromgren indices (b — y), mi and c\ were used to 
derive absolute visual magnitudes My (see Appendix B). 
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If the Hipparcos parallax (ESA 1997) of the star is avail- 
able with an error cr(7r)/-7r < 0.3, then My was also deter- 
mined directly and averaged with the photometric value. 
The bolomctric correction was adopted from Alonso et al. 
(1995) and the stellar mass derived by interpolating in 
the My- log T e ff diagram between the a-element enhanced 
evolutionary tracks of VandenBerg et al. (2000). As dis- 
cussed in Paper I, the error of log g is estimated to be about 
±0.15 dex. 

For one star, CD — 71°1234, neither uvby-f3 photometry 
nor the Hipparcos parallax is available. In this case, we 
estimated log g by requiring that the difference in the Fe 
abundance derived from Fe I and Fe II lines should equal the 
average difference between the two sets of Fe abundances 
for the other stars (see Sect. 4.1). 

The derived values of T e g and log g are mutually depen- 
dent and are also affected by [Fe/H]. Hence, the procedure 
of determining atmospheric parameters and Fe abundances 
was iterated until consistency was achieved. The final val- 
ues obtained are given in Table 1. 



4. Stellar abundances 

Using MARCS model atmospheres with the parameters 
listed in Table 1, the Uppsala eqwidth and bsyn programs 
were used to compute equivalent widths of observed lines 
as a function of the corresponding element abundance. By 
interpolation to the observed equivalent width, the abun- 
dance of the element is then derived. A basic assumption 
is LTE, but the effect of possible departures from LTE will 
be discussed. 



4.1. Iron 

The Fe abundances are based on the Fen lines listed in 
Table 2 of Paper I. The gf values of these lines were derived 
by an inverted abundance analysis of four 'standard' stars, 
for which iron abundances were adopted from Nissen et al. 
(2002). Hence, our [Fe/H] scale is tied to that of Nissen et 
al., who used a set of very weak Fen lines in the red spec- 
tral region to determine differential Fe abundances with re- 
spect to the Sun. The gf values agree very well with those 
recently derived by Melendez et al. (2006) from Fe n radia- 
tive lifetime measurements of Schnabel et al. (2004) and 
branching ratios computed by R.L. Kurucz. The average 
difference of log gf (Melendez et al. — Paper I) for the 19 
Fen lines is +0.04 dex with a scatter of ±0.06 dex only. 

Collisional line broadening was included in accordance 
with Barklcm & Aspelund- Johansson (2005). For the large 
majority of stars, the Fen lines are so weak that the de- 
rived metallicity is insensitive to possible errors in the 
van der Waals damping constant and also practically in- 
dependent of the microturbulence parameter. Stars with 
[Fe/H] < — 1.5 are assumed to have £turt>= 1.5kms _1 . For 
more metal-rich stars £turb was determined by requesting 
that the derived [Fe/H] value should be independent of 
equivalent width. 

For all new stars, the Fe abundance was also deter- 
mined from seven unblended Fei lines selected to be of 
similar strengths as the Fen lines. Wavelengths, (//-values 
(adopted from O'Brian et al. 1991) and measured equiva- 
lent widths are given in Table A.l. Collisional broadening 
data were adopted from Barklem et al. (2000b). 



0.3 
0.2 

,— , 0.1 
X 

0.0 



-0.1 
-0.2 
-0.3 



^ " * J, " •* 



5600 



5800 



6000 6200 
T M (K) 



6400 



6600 



<i 



0.3 
0.2 
0.1 
0.0 
-0.1 
-0.2 
-0.3 



-3.0 



-2.5 



-2.0 
[Fe/H] 



1.5 



-1.0 



Fig. 5. The difference in S abundances derived from the A9212.9 
and A9237.5 Si lines as a function of T e g and [Fe/H]. The rms 
deviation between the two sets of abundances is ±0.06 dex. 



The average difference between Fe abundances derived 
from Fell and Fei lines is 0.14 dex with a rms deviation of 
±0.03 dex. This small scatter testifies to the high internal 
precision obtained. The systematically higher Fe abundance 
obtained from Fe n lines may be due to non-LTE effects on 
the Fei lines (see discussion in Asplund 2005), although 
the scale of T e g and/or the gf- values could also play a role. 
We note that the Fe n - Fe I differences are very similar to 
those estimated in Asplund et al. (2006). The [Fe/H] values 
given in Table 1 are based on the Fe II lines, and refer to an 
adopted solar iron abundance of loge(Fe) Q = 7.50, which is 
close to the value inferred from weak Fe II lines in the solar 
spectrum using a MARCS model for the Sun (Nissen et al. 
2002). 

4.2. Sulphur 

Sulphur abundances were derived from the observed equiv- 
alent widths of the high excitation Si lines at 8694.6, 
9212.9 and 9237.5 A. In two of the most metal-rich stars 
(HD 179626 and HD 194598), the A8694.0 S i line could also 
be detected, and hence its corresponding S abundance was 
averaged with that from the A8694.6 line. The atomic data 
for the lines are given in Table A.l. The gf- values (adopted 
from Coulomb approximation calculations of Lambert & 
Luck 1978) agree well with the experimental values of 
Bridges & Wiese (1967). Somewhat different (//-values have 
been calculated by Biemont et al. (1993) but, as noted in 
Paper I, our adopted g /-values lead to a solar abundance 
loge(S) = 7.20, when the AA8694.0, 8694.6 lines are an- 
alyzed with a MARCS model of the Sun. This is close to 
the meteoritic abundance of sulphur, loge(S) = 7.16 ±0.04 
(Asplund et al. 2005). 

As mentioned in Sect 2, the equivalent widths of the S I 
lines at 9212.9 and 9237.5 A were measured relative to the 
flux in the wings of the Paschen-zeta H I line at 9229A. In 
Paper I, the effect of P( was neglected, but here we include 
the opacity contribution from the wings when computing 
the equivalent widths. As in the case of H/3, a version of the 
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Fig. 6. UVES observations (dots) of the AA9212.9, 9237.5 S I 
lines in the spectrum of HD 181743. The four synthetic profiles 
marked by a, b, c and d refer to loge(S)= 5.33, 5.63, 5.93 and 
6.23 corresponding to [S/Fe] = 0.0, 0.3, 0.6 and 0.9, respectively. 



Uppsala bsyn program kindly supplied by P. Barklem was 
used in these calculations. The maximum effect of P£ on the 
derived S abundances occurs for the hottest turnoff stars 
and amounts to an increase of loge(S) by about 0.07 dex 
in the case of the 9237.5 A line and about 0.02 dex for the 
9212.9 A line. After inclusion of the P£ opacity, there is an 
excellent agreement between S abundances derived from 
the two lines with no significant trend of the deviation as a 
function of T e g or [Fe/H] (see Fig. 5), in contrast to the cor- 
responding figure in Paper I, where a small but significant 
trend with T e g was seen. 

The sulphur abundances derived from the equivalent 
widths of the three S I lines are given in Table 1 together 
with the average value of [S /Fe] . For some stars it was tested 
that a detailed synthesis of the observed profiles of the S I 
lines yields practically the same abundances as the equiva- 
lent widths (see example in Fig. 6). 

Table 1 includes three bright halo stars, HD 84937, 
HD 140283 and HD 181743, claimed by others to have very 
high sulphur abundances, [S/Fe] £ + 0.6, whereas we find 
these stars to have [S/Fe] around +0.3 dex. In trying to 
find the reason for this discrepancy, we briefly discuss these 
stars. Fig. 7 shows the spectrum synthesis of their S I lines 
at 8694 A. The instrumental and stellar line broadening 
profile was approximated by a Gaussian with a FWHM 
of Skins" 1 in the case of HD 84937 and 6kms -1 for the 
other two stars. These values were derived from fitting the 
Fei line at 8688.6 A. 

4.2.1. HD84937 

From a KECK/HIRES spectrum of HD 84937, Takada- 
Hidai et al. (2002) measured the equivalent width of the 
A8694.6 Si line to be W x = 3.4mA and derived loge(S)= 
5.7 corresponding to [S/Fe]= +0.6. From our UVES spec- 
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Fig. 7. UVES spectra (jagged line) of the region around the 
A8694.6 S I line. For each star, the three synthetic profiles 
marked by a, b and c correspond to [S/Fe] = 0.3, 0.6 and 0.9. 



trum we measure W\ = 1.6 ± 0.3 mA, where the quoted (1- 
sigma) error is estimated from the S /N of the spectrum and 
the uncertainty of the continuum setting. The correspond- 
ingly derived sulphur abundance is loge(S) = 5.40 ± 0.08. 
As shown in Fig. 7, this agrees well with the spectrum syn- 
thesis of the line; a value as high as loge(S)= 5.7 is clearly 
excluded. Furthermore, we obtain loge(S)= 5.39 and 5.48 
from the A9212. 9, 9237.5 Si lines. Altogether, our sulphur 
abundance of HD 84937 is about a factor of two lower than 
the value found by Takada-Hidai et al. (2002). 

4.2.2. HD 140283 

Using the UVES Paranal Observatory Project (Bagnulo et 
al. 2003) spectrum of HD 140283, Takeda et al. (2005) de- 
rived [S/Fe] ~ +1.0 (corresponding to loge(S)^ 5.8) from 
the A8694.6 Si line (see Fig. 7 in Takeda et al. 2005). As 
seen from our Fig. 7, we barely detect this line in our S/N 
- 600 spectrum of HD 140283; clearly, loge(S) is much 
lower than 5.8 dex. From the measured equivalent width 
of the A8694.6 line, W x = 0.7 ± 0.3 mA, we get loge(S) 
= 5.17 ± 0.20 in good agreement with the result from the 
AA9212.9, 9237.5 Si lines, i.e. loge(S) = 5.07 and 5.06, re- 
spectively. Judging from Fig. 7 in Takeda et al. (2005), the 
UVES POP spectrum of HD 140283 is affected by residual 
fringing, which in this case apparently has led to a spu- 
rious enhancement of the A8694.6 S I line, as also pointed 
out by Korn & Ryde (2005). Furthermore, it seems that 
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Takeda et al. have placed the continuum in the A8694.6 re- 
gion too high. Using the same UVES POP spectrum with 
a more reasonable continuum, Takada-Hidai (2005) derived 
loge(S) = 5.53 from the A8694.6 line. This is, however, still 
much higher than our value and also higher than the value 
of loge(S) = 5.18, which they derive from the A9237.5 Si 
line. 



4.2.3. HD 181743 

For this star, Caffau et al. (2005) obtained loge(S)= 6.23 
and [S/Fe] = +0.84. The result is based on an apparently 
clear detection of the A8694.6 S I line in a UVES spectrum 
with S/N = 180 (P. Bonifacio, private communication). 
However, we cannot detect this line in our S/N ~ 280 spec- 
trum as seen from Fig. 7, but estimate a 3-sigma upper 
limit of the equivalent width, W\ < 2 mA corresponding 
to loge(S) < 5.8. Caffau et al. did not use other Si lines 
for this star, whereas we derive loge(S)= 5.62 from the 
AA9212.9, 9237.5 lines. As seen from Fig. 6 the fit of the 
corresponding synthetic profiles to the observed lines is ex- 
cellent. 

We note that the large differences between sulphur 
abundances derived in the cited works and in the present 
paper cannot be explained in terms of different values of 
T c ff and logg. In our view, it is more likely that the S 
abundances derived by others from the weak A8694.6 Si 
line have been overestimated due to fringing residuals in 
their spectra. In general, more precise S abundances can 
be derived from the AA9212.9, 9237.5 lines as they are more 
easily detected and less susceptible to observational uncer- 
tainties. In this connection, it should be noted that Caffau 
et al. find a few stars to have very high [S/Fe] values even 
on the basis of these lines. Most remarkable is BD +02°263, 
for which they get [S/Fe] = +0.91. This star is, however, 
a single-lined spectroscopic binary according to Latham et 
al. (2002), which may have affected the determination of 
T e g and log g and hence the abundances of S and Fe. 

4.2.4. 3D and non-LTE effects on the Si lines 

In Paper I it was shown that the application of 3D hydrody- 
namical model atmospheres instead of classical ID models 
has a rather small effect on the derived sulphur abundance 
(see also Asplund 2005). For the large majority of our stars, 
3D models lead to an increase of loge(S) by 0.05 to O.lOdex. 
This increase is, however, compensated by about the same 
increase in iron abundance (see Table 5 in Paper I for de- 
tails). Hence, [S/Fe] is practically unchanged. 

It was noted in Paper I that the good agreement be- 
tween S abundances derived from the weak A8694.6 S I line 
and the stronger AA9212.9, 9237.5 lines suggests that non- 
LTE effects on the sulphur lines are small. Thanks to the 
recent non-LTE calculations of Takeda et al. (2005), this 
statement may now be quantified. From extensive statis- 
tical equilibrium calculations Takeda et al. show that the 
AA9212.9, 9237.5 lines can suffer from significant negative 
non-LTE corrections, whereas the non-LTE effects on the 
A8694.6 line are small. The size of the corrections depend 
on the rate of inelastic collisions with neutral hydrogen, for 
which Takeda et al. adopt the classical approximation of 
Drawin (1968, 1969) in the version of Steenbock & Holweger 
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Fig. 8. Comparison of S abundances derived from the S I 9212.9, 
9237.5 A pair and the A8694.6 line. Filled circles refer to LTE. 
Open circles include non-LTE corrections according to Takeda 
et al. (2005) with the hydrogen collisional parameter Sh = 1, 
and crosses refer to Sh = 0.1. 



(1984) with a scaling factor Sh that is varied from 10 3 to 
10. 

As seen from the tables of Takeda et al. (2005), the non- 
LTE effects vary quite strongly with T e g, log g and [Fe/H], 
but interpolation in the tables allows us to derive non-LTE 
corrections for our set of stars. In Fig. 8, the resulting differ- 
ence of the S abundance derived from the A9212.9, 9237.5 A 
pair and the A8694.6 line has been plotted for two values of 
Sh and also in the LTE case. As seen, both LTE and Sh = 1 
give a satisfactory agreement between the two sets of S 
abundances, whereas the case of Sh = 0.1 can be excluded. 
The average difference and dispersion of log e(S)g2i2, 9237 ~ 
loge(S) 8694 are as follows: +0.04 + 0.10 (LTE), -0.06 + 0.09 
(non-LTE with S H = 1), and -0.27 + 0.08 (non-LTE with 
Sh = 0.1). This suggests that hydrogen collisions are quite 
efficient in thermalizing the S 1 atoms. Still, the non-LTE 
effects for the AA9212.9, 9237.5 lines may be of importance 
for the derived [S/Fe] trend, due to varying non-LTE abun- 
dance corrections as a function of [Fe/H]. In the case of 
Sh = 1, the corrections obtained by Takeda et al. range 
from about — 0.06dex for the coolest of our stars to about 
— 0.25dex for the hottest and most metal-poor stars (e.g. 
G 64-12). 

4.2.5. CRIRES observations of the 1.046 Si triplet 

As part of the science verification of the ESO/VLT 
Cryogenic high-resolution IR Echelle Spectrograph, 
CRIRES (Kauri et al. 2004), a spectrum around the 
1.046 /j,m Si triplet was obtained for G 29-23 (V = 10.19, 
[Fe/H]= -1.69) on October 6, 2006. The entrance slit 
width of CRIRES was set at 0.4 arcsec, which corresponds 
to a resolution of A/ A A ~ 50 000 with four detector 
pixels per spectral resolution bin AA. In order to improve 
removal of sky emission and detector dark current, the 
observations were performed in nodding mode with a shift 
of 10 arcsec between the two settings of the star in the 
slit. The exposure time was 2400 sec. The seeing was 
rather poor (about 1.3 arcsec) but adaptive optics was 
applied to improve the stellar image, and the combined 



P.E. Nissen et al.: Sulphur and zinc abundances in Galactic halo stars revisited 



9 



Table 2. Atomic data for the 1.046 fj,m Si triplet. Column 5 
lists the equivalent width in the spectrum of G 29-23 and the 
last column gives the derived LTE sulphur abundance. 





Wavelength 

A 


Exc. Pot. 
eV 


log gf 


W A (G29-23) 
mA 


loge(S) 


Si 


10455.5 


6.86 


0.26 


39.5 


5.82 


Si 


10456.8 


6.86 


-0.44 


10.1 


5.75 


Si 


10459.4 


6.86 


0.03 


25.7 


5.77 



spectrum has a very satisfactory S/N of 330 per spectral 
dispersion pixel. For comparison, we note that our UVES 
spectrum of G 29-23 (exposure time 1800 sec) has S/N 
~ 200 around the Si lines at 9212-9238 A. 2 Furthermore, 
unlike the UVES near-IR spectrum, the CRIRES spectrum 
at 1.046 /im is not plagued by telluric lines and fringing 
residuals. 

The CRIRES data was reduced by using standard IRAF 
tasks for subtraction of bias and dark, spectrum extraction, 
and wavelength calibration. Note that CRIRES supports 
high accuracy wavelength calibration by means of a Th- 
Ar hollow cathode lamp (Kerber et al. 2006). Flatfielding 
was performed by the aid of the extracted spectrum of the 
bright (V = 4.28) B9III star HD 15315. The resulting spec- 
trum of G 29-23 is shown in Fig. 9 and compared with a 
spectral synthesis of the three sulphur lines. 

Atomic data for the 1.046 /im Si triplet are given in 
Table 2 together with the measured equivalent widths in 
the G 29-23 spectrum. As in the case of the other S I lines 
the gf values were taken from Lambert & Luck (1978) 
and collisional broadening data are from Barklem et al. 
(2000b). The derived LTE S abundances given in Table 2 
agree very well with the S abundances derived from the S I 
lines observed with UVES. An average of [S/Fe] = +0.29 
is obtained from the UVES lines and [S/Fc] = +0.27 from 
the CRIRES lines. For all six Si lines the mean sulphur 
abundances is loge(S)= 5.79 and the rms deviation is only 
0.035 dex. Although one should not put too much weight 
on a single star, we consider this good agreement to be an 
important check of the reliability of our sulphur abundance 
determinations . 

According to Takeda et al. (2005), the non-LTE cor- 
rections for the 1.046 /jm Si triplet are somewhat smaller 
than in the case of the A9212.9, 9237.5 A lines. Hence, one 
might in principle also use the 1.046 /im triplet to study 
non-LTE effects and calibrate the S*h parameter like we did 
when comparing S abundances from the A8694.6 line and 
the AA9212.9, 9237.5 pair. This would, however, require IR 
observations for many more stars to provide statistically 
significant results. 

4.3. Zinc 

The Zn abundances are derived from the equivalent widths 
of the two Zn I lines at 4722.2 and 4810.5 A. Adopted gf val- 
ues (see Table A.l) are from Biemont & Godefroid (1980), 
and collisional broadening data were taken from Barklem et 
al. (2000b). In Paper I, a solar Zn abundance of loge(Zn) Q = 

2 The UVES spectra also cover the 1.046 /im Si triplet but 
due to the low efficiency of the MIT CCD at this wavelength, 
the quality is far too low (S/N < 25) for a precise determination 
of the sulphur abundance. 
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Fig. 9. The CRIRES spectrum of G 29-23 around the 1.046 fim 
S I triplet (dots) compared with synthetic profiles for three sul- 
phur abundances, loge(S) = 5.51, 5.81 and 6.11 corresponding 
to [S/Fe] = 0.0, 0.3 and 0.6, respectively. As seen, the [S/Fe] = 
0.3 case provides an excellent fit to the observations. 

4.57 was derived from the two lines, but this value is sensi- 
tive to the solar microturbulcncc parameter. Here, wc adopt 
loge(Zn) = 4.60 as derived by Biemont & Godefroid (1980) 
from six Zni lines. This value is in good agreement with 
the meteoritic value, loge(Zn) = 4.61 ± 0.04 (Asplund et 
al. 2005). 

As seen from Table A.l, the zinc lines are very weak 
in turnoff stars with metallicities below —2. The equiv- 
alent widths are on the order of 2-3 mA in stars with 
[Fe/H] ~ —2.5. Still, the Zn lines can be clearly detected at 
metallicities below —3 as seen from Fig. 10. This is con- 
nected to the fact that the Zn/Fe ratio appears to rise 
steeply below [Fe/H] ~ — 3. 

Takeda et al. (2005) have also performed non-LTE cal- 
culations for the AA4722. 2, 4810.5 Zni lines. Generally, the 
absolute value of the non-LTE corrections are smaller than 
in the case of the AA9212.9, 9237.5 Si lines, and they have 
opposite sign. For a typical metal-poor turnoff star, the 
non-LTE correction of the zinc abundance is on the or- 
der of +0.10 dex, whereas the correction is —0.25 dex for 
the sulphur abundance. For cooler and more metal-rich 
stars the non-LTE corrections on zinc are less than 0.1 dex. 
These corrections refer to a hydrogen collisional parameter 
Sh = 1, and increase somewhat if Sh < 1. 

As discussed in Paper I, the absolute values of the 3D 
corrections of the Zn abundances are less than 0.1 dex and 
they are similar to the 3D corrections of the Fe abundances 
based on Fen lines. Hence, the derived Zn/Fe ratio would 
not change significantly (at least not under the assumption 
of LTE) if 3D models were applied. 

5. Results and discussion 

5.1. S/Fe vs. Fe/H 

The derived sulphur abundances arc given in Table 1 . In the 
calculation of [S/Fe], solar abundances loge(S) Q = 7.20 and 
loge(Fe) Q = 7.50 were adopted. Had we instead applied the 
newest solar abundances, loge(S)o= 7.14 and loge(Fe)© = 
7.45, determined by Asplund et al. (2005), the [S/Fe] values 
would have increased by only 0.01 dex. 
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Fig. 10. The A4810.5 Znl line in UVES spectra of three stars. 
The spectra of G 04-37 and CD —24° 17504 have been offset in 
relative flux by 0.04 and 0.08, respectively. 



The [S/Fc] values in Table 1 are plotted vs. [Fe/H] in 
Fig. 11 together with data for 25 disk stars from Chen et 
al. (2002). The error bars shown refer to the 1-a statistical 
error caused by the uncertainty of the observed equivalent 
widths and the atmospheric parameters of the stars. In the 
case of [Fe/H], which is based on many Fen lines, the dom- 
inating error comes from the uncertainty of the gravity; 
cr(log.g) = 0.15 dex corresponds to cr[Fe/H] = 0.05 dex (see 
Table 4 in Paper I). The log g induced error of [S/Fe] is, 
on the other hand, negligible, because a change in grav- 
ity has nearly the same effect on S and Fe abundances. 
In the case of [S/Fe] the major error contribution comes 
from the uncertainty of T e g and the error of the equivalent 
width measurements. For the three most metal-poor stars 
with [Fe/H] < —3, the W\ induced error becomes particu- 
lar large, because only the A9212.9 S I line is strong enough 
to be detected and it is as weak as W\ ~ 5 mA. 

As seen from Fig. 11, the LTE values of [S/Fe] suggest 
a small slope of [S/Fe] as a function of [Fe/H], whereas the 
non-LTE values give a nearly flat relation. The mean [S/Fe] 
and the rms deviation for the halo stars are: 

< [S/Fe] > = 0.319 ±0.078 in the LTE case, and 

< [S/Fe] > = 0.208 ±0.067 in the non-LTE case. 

We note that the four single- lined spectroscopic binaries 
(see Table 1) do not show any particular large or systematic 
deviations from the mean value of [S/Fe]. 

The trend of [S/Fe] in Fig. 11 corresponds rather well 
to analogous trends recently derived for [Mg/Fe], [Si/Fe] 
and [Ca/Fe] (Jonsell et al. 2005, Gehren et al. 2006). 
Qualitatively, our data also agree with the trend of [S/Fe] 
predicted from models of the chemical evolution of our 
Galaxy. Based on the Woosley & Weaver (1995) yields 
of Type II SNe, Goswami & Prantzos (2000) obtain a 
near-constant level [S/Fe] ~ ±0.4 in the metallicity range 
-3 < [Fc/H] < -1 and a decline of [S/Fe] for [Fe/H] > -1 
due to the additional supply of iron from Type la SNe. 
Chiappini et al. (1999), on the other hand, obtain a slope 
of [S/Fe] in the range —3 < [Fe/H] < — 1 like the one we 
determine in the LTE case, because they assume that Type 




Fig. 11. [S/Fe] vs. [Fe/H] for our sample of halo stars supple- 
mented with disk stars (shown as crosses) from Chen et al. 
(2002). Filled circles with error bars show data based on S abun- 
dances derived from the AA9212.9, 9237.5 S I lines, and open 
circles show data based on the weak A8694.6 S I line. In the 
upper panel LTE has been assumed in deriving the S abun- 
dances, whereas the lower panel includes non-LTE corrections 
from Takeda et al. (2005) with the hydrogen collisional param- 
eter Sh = 1. In both cases, the iron abundances are the LTE 
values derived from Fell lines. 



Ia SNe start contributing with iron already in the Galactic 
halo phase. Finally, Kobayashi et al. (2006) predict a level 
of [S/Fe] ~ ±0.4 based on new yield calculations in a model 
where hypernovae are assumed to have the same frequency 
as Type II SNe. We note that our data do not allow us 
to establish conclusively whether a small slope of [S/Fe] is 
present or not due to the uncertainty about the non-LTE 
corrections. It should also be noted that the predicted level 
of [S/Fe] among halo stars depends critically on the as- 
sumed mass cut between the ejecta and the collapsing core 
in massive SNe. Hence, any detailed comparison between 
observed and predicted trends of [S/Fe] seems somewhat 
premature. 

It is interesting to compare the observed scatter of 
[S/Fe] with the scatter predicted from stochastic chemical 
evolution models. In sharp contrast to the results of Caffau 
et al. (2005), we find a very small scatter (±0.07 dex) of 
[S/Fe] among the halo stars, and there is no indication 
that the scatter increases significantly towards the most 
metal-poor stars. This scatter is not much higher than ex- 
pected from the observational errors of [S/Fe]. Hence, the 
cosmic scatter in [S/Fe] at a given [Fe/H] must be less than 
0.07 dex. A similar low scatter has been obtained for other 
a-capture elements, when high precision and homogeneous 
data have become available. Nissen et al. (1994) found the 
scatter of Mg, Ca and Ti relative to Fe to be < 0.06 dex 
in the metallicity range —3 < [Fe/H] < —1.5. More re- 
cently, Arnone et al. (2005) found the cosmic scatter of 
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[Mg/Fe] to be < 0.06 dex for 25 turnoff stars in the range 
-3.2 < [Fc/H] < -2 3 . At still lower metallicities, Cayrel 
et al. (2004) find a scatter of 0.10 - 0.15 dex for [Mg/Fe], 
[Si/Fc], [Ca/Fe] and [Ti/Fe] for 30 giant stars ranging in 
[Fc/H] from —4.1 to —2.7. A similar low scatter of these 
ratios has been derived by Cohen et al. (2004) for 28 dwarf 
stars in the metallicity range —3.6 < [Fe/H] < —2.0. 

The observed scatter of a-capture elements relative to 
Fe is much smaller than predicted from stochastic models 
of the chemical evolution of metal-poor systems (Argast et 
al. 2000, 2002, Karlsson & Gustafsson 2005). This is con- 
nected to the fact that calculated yields of Type II SNe vary 
strongly with progenitor mass. According to Nomoto et al. 
(1997), the yield ratio S/Fe increases by a factor of 12 when 
the progenitor mass is changed from 13 to 40 solar masses. 
Hence, in the early Galaxy, where only a few supernovae 
enrich the interstellar gas of a star-forming cloud according 
to current models, one would expect a much higher scatter 
in [S/Fe] than the derived value of < 0.07 dex for our sam- 
ple of stars. Sulphur has not been specifically modelled in 
any of the stochastic models, but in the case of Si that has 
similar yield variations as S, Argast et al. (2000) predict 
a rms scatter of ±0.35 dex at [Fe/H] = -3.0, ±0.25 dex 
at [Fc/H] = -2.5, and ±0.12 dex at [Fe/H] = -2.0. As 
discussed in detail by Arnone et al. (2005) and Karlsson 
& Gustafsson (2005), possible explanations of this discrep- 
ancy are: i) The calculated yield ratios are wrong, i.e. the 
released amount of alpha-elements relative to Fe is nearly 
independent of progenitor mass, ii) the IMF is biased to a 
narrow mass range in the early Galaxy, and in) the mixing 
of SNe ejecta is much more rapid than assumed in the mod- 
els, such that a large number of SNe always contribute to 
the enrichment of a cloud. The last possibility means that 
the mixing time scale of supernova ejecta is considerably 
shorter than the cooling time of star-forming gas clouds. 

It should be noted that the scatter in the abundance 
of a-capture elements relative to Fe probably increases 
towards the more metal- rich end of the halo. Paper I in- 
cludes six halo stars with [Fe/H] > — 1, four of which have 
[S/Fe] ~ +0.3 and two have [S/Fe] ~ 0.0. This scatter is 
in accordance with the results of Nissen & Schuster (1997), 
who find a similar dichotomy in [a/Fe] for halo stars with 
— 1 < [Fe/H] < —0.6. As suggested in their paper, the ex- 
planation of the scatter may be that the metal-rich a-poor 
stars have formed in the outer regions of the Galaxy where 
the star formation has proceeded so slowly that iron from 
Type la SNe has been incorporated in them. Supporting 
evidence has been obtained by Gratton et al. (2003), who 
find [a/Fe] to be lower and more scattered in stars belong- 
ing to the 'accretion' component of the Galaxy than in stars 
belonging to the 'dissipative' component. 



5.2. Zn/Fe vs. Fe/H 

The LTE values of [Zn/Fe] are plotted as a function of 
[Fe/H] in the upper panel of Fig. 12. In the lower panel, non- 
LTE corrections of [Zn/Fe] from Takeda et al. (2005) have 
been included. The error bars refer to the l-cr statistical 
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This range is on our metallicity scale. The metallicities of 
Arnone et al. (2005) are based on Fel lines and a T e s scale that 
is on the average ~ 200 K lower than our scale. For 14 stars in 
common with our sample, they get [Fe/H] values that are on the 
average 0.20 dex lower than our values. 



Fig. 12. [Zn/Fe] vs. [Fe/H] for our sample of halo stars. The up- 
per panel refers to LTE Zn abundances. The lower panel includes 
non-LTE corrections from Takeda et al. (2005) corresponding to 
a hydrogen collisional parameter Sn = 1. 



error caused by the uncertainty of the observed equivalent 
widths and the atmospheric parameters of the stars. For 
the three most metal-poor stars, the HVinduced error is 
large, because the equivalent widths of the Zni lines are 
only about 1 mA. 

As seen from the upper panel of Fig. 12, the LTE val- 
ues of [Zn/Fe] are close to zero in the metallicity range 
-2 < [Fc/H] < -1. In the range -2.7 < [Fc/H]' < -2, 
all [Zn/Fe] values are, however, positive with an average 
value of [Zn/Fe] ~ +0.1. Furthermore, [Zn/Fe] seems to rise 
steeply below [Fe/H] = — 3 to a value of [Zn/Fe] ~ +0.5. 
This trend is reinforced if non-LTE corrections correspond- 
ing to Sh = 1 arc applied, because the corrections to 
[Zn/Fe] are positive and increase with decreasing [Fe/H]. 
Unlike the case of sulphur, we had no possibility to calibrate 
the Sh parameter for zinc (our two Zni lines have nearly 
the same excitation potential), but we note that if Sh < 1 
(corresponding to a lower efficiency of thermalizing the Zn 
atoms by inelastic collisions with hydrogen atoms) then the 
non-LTE corrections for the metal-poor stars would be still 
larger. 

The trend of [Zn/Fe] in Fig. 12 (LTE case) corresponds 
very well to the trend found by Cayrel et al. (2004) from 
an LTE analysis of the AA4722.2, 4810.5 Zni lines in UVES 
spectra of 35 giant stars with —4.2 < [Fe/H] < —2. Hence, 
there can be little doubt that Zn is indeed overabundant 
with respect to iron for very metal-poor stars, especially 
because the possible non-LTE corrections go in the direc- 
tion of making the overabundance larger. 

Zinc is a key element in recent studies of nucleosyn- 
thesis and chemical evolution in the early Galaxy (Umcda 
& Nomoto 2003, 2005; Nomoto et al. 2006; Kobayashi et 
al. 2006). As shown by Kobayashi et al. (2006), traditional 
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yields of Type II SNe (Nomoto et al. 1997) correspond to 
[Zn/Fc] ~ —1, i.e. far below the observed values of [Zn/Fe]. 
In order to explain a level of [Zn/Fe] around zero, one has 
to invoke models of hypernovae, i.e. core collapse SNe with 
explosion energy E ^ 10 52 ergs. Furthermore, a mixing and 
fallback mechanism or an asymmetric explosion (Maeda 
& Nomoto 2003) has to be introduced in order to bring 
up sufficient Zn from layers with complete Si burning. By 
including such hypernovae with a frequency of 50% rela- 
tive to Type II SNe and a Salpeter IMF, Kobayashi et al. 
(2006) predict [Zn/Fe] ~ 0.1 in the Galactic halo. The up- 
turn of [Zn/Fe] at the lowest metallicities is, however, not 
predicted. In order to explain [Zn/Fe] ^ + 0.3, Nomoto et 
al. (2006) suggest that stars with [Fe/H] < —3 have been 
formed from the ejecta of Pop. Ill hypernovae with very 
large explosion energy. Another possibility for high Zn/Fe 
production is from core-collapsing, very massive stars with 
M ~ 500 - 1OOOM (Ohkubo et al. 2006). 

5.3. S/Zn vs. Zn/H 

As explained in the Introduction, the S/Zn ratio is poten- 
tially of great value in the analysis of the metallicities of 
distant galaxies detected via the absorption lines they pro- 
duce in the spectra of quasars; the damped Lyman-alpha 
systems (DLAs) are a particularly important subset of such 
absorbers (Wolfe at al. 2005) . The role of S and Zn in these 
studies stems from the fact that both elements have a low 
condensation temperature and are not expected to be sig- 
nificantly depleted onto dust. This is not the case, for ex- 
ample, for refractory elements such as Si and Fe. Indeed, 
one of the motivations for the present study was to use the 
behaviour of [S/Zn] vs. [Zn/H] in Galactic stars as a guide 
to the interpretation of the star formation history of high 
redshift galaxies. 

Figure 13 shows [S/Zn] vs. [Zn/H] for our halo stars 
except the three with [Fe/H] < — 3, which have such large 
errors in both S and Zn that the ratio S/Zn becomes very 
uncertain. In addition, we have included six halo stars with 
— 1.1 < [Zn/H] < —0.6 from Paper I and 14 disk stars with 
S and Zn abundances from Chen et al. (2002; 2004). The 
upper panel shows LTE values; in the lower panel non-LTE 
corrections (Sh = 1) from Takeda et al. (2005) have been 
applied. By comparing the two panels, the importance of 
the non-LTE corrections is readily apparent. When the cor- 
rections are applied, there is an overall decrease in [S/Zn] 
at all but the highest values of [Zn/H] considered here. 
Furthermore, non-LTE effects are most significant at low 
metallicities with the result that, apparently, [S/Zn] reverts 
to solar values when [Zn/H] ^ — 2. Such behaviour is un- 
usual but, given our current limited understanding of the 
nucleosynthesis of Zn, cannot be excluded. 

Taken at face value, the lack of a strong metallicity trend 
in the lower panel of Fig. 13 would indicate that the use- 
fulness of the S/Zn ratio as a 'clock' of the star- formation 
history is rather limited — departures from the solar value 
are generally small and, in particular, not much greater 
than the typical measurement error in DLAs. In Table C.l 
we have collected all published, or about to be published, 
measurements of [S/Zn] in DLAs secured with echelle spec- 
trographs on 8-10 m class telescopes. In comparison with 
the corresponding table in Paper I, the number of DLAs 
with accurate measures of the abundances of sulphur and 
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Fig. 13. [S/Zn] vs. [Zn/H] for Galactic stars and damped 
Lyman-alpha systems. Typical 1-a error bars are shown. 



zinc has doubled since 2004, from 10 to 20. These data are 
shown as triangles in Fig. 13. 

Whereas in previous papers (e.g. Paper I) a system- 
atic difference had been noted between the [S/Zn] ratio 
in Galactic metal-poor stars and in DLAs of comparable 
metallicities (evident in the top panel of Fig. 13), the sit- 
uation is considerably less clear-cut if one adopts the non- 
LTE corrections to the abundances of S and Zn. If these 
corrections are appropriate, it would explain why so few 
DLAs show enhanced [S/Zn] values, when other ratios of 
alpha-capture to iron-peak elements apparently do, once 
dust depletions are accounted for (e.g. Prochaska & Wolfe 
2002). A minor difference - which needs to be quantified 
with larger datasets - between Galactic stars and DLAs is 
the larger scatter in the [S/Zn] ratio possibly exhibited by 
the latter. It is intriguing, in particular, to find an under- 
abundance of S relative to Zn in some DLAs with relatively 
high overall metallicities. If confirmed, the larger disper- 
sion of the ratio in DLAs could be an indication that, as 
one might expect, the host galaxies of DLAs experienced a 
variety of different star formation histories, and don't nec- 
essarily follow closely the abundance trends seen in a single 
galaxy like the Milky Way. 



6. Summary and conclusions 

High S/N UVES spectra of 40 metal-poor halo stars have 
been used to derive S, Fe and Zn abundances. For one 
star we also present novel observations of the Si triplet 
at 1.046 ^m carried out with the ESO VLT CRIRES spec- 
trograph. These data confirm the S abundance obtained 
from other S I lines and demonstrates the high efficiency of 
CRIRES in obtaining high-resolution infrared spectra. 

From the results presented in this paper, we conclude 
that [S/Fe] in Galactic halo stars shows the same kind of 
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dependence on [Fe/H] as [Mg/Fe], [Si/Fe] and [Ca/Fe], i.e. 
a near-constant ratio at a level of +0.2 to +0.3 dex. This 
strongly suggests that sulphur in the Galactic halo was 
made by a-capture processes in massive SNe, as predicted 
from current models of Galactic chemical evolution with 
yields calculated for hydrostatic and explosive oxygen and 
silicon burning. 

Among the 40 halo stars observed, we do not find a sin- 
gle case of the high S/Fe ratios ([S/Fe] > 0.60) claimed in 
some recent papers (Israelian & Rebolo 2001;Takada-Hidai 
et al. 2002; Caffau et al. 2005). The majority of these high 
values were based on the very weak A8694.6 S I line and, 
as shown in Sect. 4.2, our very high S/N, fringe-corrected 
spectra do not confirm the claimed strength of this line 
in three metal-poor main-sequence stars. A more convinc- 
ing detection of the A8694.6 S I line has been obtained by 
Israelian & Rebolo (2001) for two giant stars, HD 2665 and 
HD2796 with [Fe/H] = -2.0 and -2.3. The derived [S/Fe] 
values are +0.69 and +0.81 dex, respectively. Further stud- 
ies of S abundances in such metal-poor giant stars would 
be interesting, especially if one could use the forbidden sul- 
phur line at 1.082 /zm, which is expected to be insensitive 
to non-LTE effects. 

From an observational point of view, the 
AA9212.9, 9237.5 pair of Si lines or the Si triplet at 
1.046 itm provide more precise values of the S abundance 
than the weak A8694.6 Si line, but according to Takeda 
et al. (2005) these stronger lines are quite sensitive to 
departures from LTE. The size of the non-LTE effects 
depends on the efficiency of inelastic collisions with neutral 
hydrogen. Adopting the classical formula of Drawin (1968, 
1969) with a scaling parameter Su, a comparison of S 
abundances from the A8694.6 line and the AA9212.9, 9237.5 
pair points to Su £ 1, which means that we are not too 
far away from LTE. Nevertheless, the non-LTE corrections 
may be of importance for the derived trend of [S/Fe] as 
seen from Fig. 1 1 . Hence, it would be important to perform 
improved quantum mechanical calculations for inelastic S 
+ H collisions as has been achieved for Li + H collisions 
(Belyaev & Barklem 2003). 

An interesting and robust result of our investigation is 
the small scatter of [S /Fe] at a given metallicity. It is found 
to be around ±0.07 dex in the whole range —3.2 < [Fe/H] < 
— 1. A similarly low scatter has been found for the other a- 
capture elements, i.e. [Mg/Fe], [Si/Fe], and [Ca/Fe] (Cayrel 
et al. 2004; Cohen et al. 2004). Current stochastic models 
of the chemical evolution of metal-poor systems (Argast et 
al. 2000, 2002, Karlsson & Gustafsson 2005) predict a much 
higher scatter of [a/Fe] for metallicities below —2. Clearly 
the models and/or calculated yield ratios of massive SNe 
need to be revised. 

Zinc is found to be slightly overabundant with re- 
spect to Fe in the metallicity range —3 < [Fe/H] < —2. 
Below [Fe/H] ~ —3, [Zn/Fe] increases rapidly to a level 
of +0.5 dex. Hence, our study of main-sequence halo stars 
confirms the upturn of Zn/Fe previously found for very 
metal poor giants (Primas et al. 2000; Johnson & Boltc 
2001; Cayrel et al. 2004). The high value of Zn/Fe at the 
lowest metallicities may be a signature that stars with 
[Fe/H] < —3 have been made from the ejecta of hypernovae 
with very large explosion energy as suggested by Nomoto 
et al. (2006). However, our understanding of the nucleosyn- 
thesis of Zn at all metallicities is still very limited. 



Finally, we find a much reduced 'signal' in the [S/Zn] 
ratio as a tracer of the previous history of star formation, 
if non-LTE effects to the abundances of these two elements 
are taken into account. It would be very desirable to per- 
form a check of the validity of the non-LTE corrections to 
the Zn abundance by comparing lines arising from levels 
with different excitation potentials, as we have been able 
to do here for S. However, if the non-LTE corrections ap- 
plied arc appropriate, it would explain why in general DLAs 
do not exhibit markedly supersolar [S/Zn] ratios, even at 
metallicities comparable to those of Galactic halo stars. 
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Appendix A: Line list and equivalent widths 

Table A.l gives a list of spectral lines used in this paper and 
the equivalent widths measured in the 2004 UVES spectra. 
Where no value is given, it is either because the line is too 
weak to provide a reliable abundance, or, in the case of the 
AA9212. 9, 9237.5 Si lines, is affected by residuals from the 
removal of strong telluric H 2 lines. We also note that the 
Si line at 9228.1 A is not included, because it falls close to 
the center of the Paschen-zeta line. 



Appendix B: Photometric and spectroscopic 
indices 

In Table B.l the magnitudes V and K s for our sample of 
stars are given together with the Stromgren indices (b—y), 
mi, a and the photometric index of the H/3 line, /3(phot). 
The table also gives the value of the /3(UVES) index mea- 
sured from the profile of H/3 in the UVES spectra. This 
index is defined as 

/3(UVES) = 10- p Cb +p Cr , (B.l) 

where Feb, i*br are the fluxes in two 8 A wide pseudo- 
continuum bands centered ±20 A from the H/3 line center 
and FLb, F^ r the fluxes in two 10 A wide line bands cen- 
tered ±7 A from the center of H/3. As seen from Fig. 3, 
the fluxes in the line bands are more affected by variations 
in effective temperature than the fluxes in the 'continuum' 
bands. Hence, the index becomes sensitive to T e g. 
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Table A.l. List of spectral lines and equivalent widths measured in the UVES 2004 spectra of the following twelve stars: (1) 
CD —24° 17504, (2) CD —71° 1234, (3) CS 22943-095, (4) G 04-37, (5) G 48-29, (6) G 59-27, (7) G 126-52, (8) G 166-54, (9) HD 84937, 
(10) HD 338529, (11) LP 635-14, (12) LP 651-4. 
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In Fig. B.l some of the calculated /3(UVES) values are 
plotted as a function of T e g together with the interpola- 
tion lines. As seen, the variation of /3(UVES) with \ogg 
and [Fe/H] is small compared to the change of /3(UVES) 
as a function of T c tf. For an observational error ±0.015 of 
/3(UVES) and errors ±0.15 dcx of log 5 and [Fe/H], T cS is 
determined to a precision of about ±20 K. 

As seen from Fig. B.2 there is a good correlation be- 
tween /3(phot) and /3(UVES). A maximum likelihood fit 
that takes into account the estimated errors of the two in- 
dices gives the following relation: 

/3(phot) = 1.474 + 0.129 • /3(UVES) (B.2) 

with a reduced x 2 close to one. The dominant contribution 
to the scatter in Fig. B.2 comes from the error of /3(phot). 
According to Eq. (B.2), the error of /3(UVES) (±0.015) cor- 
responds to an error of ±0.002 in /3(phot), which is a fac- 
tor of three lower than the actual error of the /3(phot) ob- 
servations (Schuster & Nissen 1988). Hence, we have used 
Eq. (B.2) to predict an improved /3(phot) value of our stars 
before calculating the interstellar reddening excess from the 
(b — y) calibration equations of Schuster & Nissen (1989) 
with a zero point shift of 0.005 mag. added (Nissen 1994). 
The so derived value of E{b — y) is given in column 10 of 
Table B.l. 

From this discussion it also follows that the error of 
T c ff determined from the photometric H/3 index would be 
three times higher than the error of T g- determined from 
/3(UVES). On the other hand, it is seen from Fig. B.3 that 
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Fig. B.l. The relation between /3(UVES) and T e s as calculated 
for a set of model atmospheres. The full drawn line corresponds 
to logp = 4.3, [Fe/H] = -2.0, the dotted line to logg = 3.9, 
[Fe/H] = -2.0, and the dashed line to log# = 4.3, [Fe/H] = 
-3.2. 



T e ff was derived from the observed value of /3(UVES) 
by quadratic interpolation between theoretical /?(UVES) 
values calculated for a grid of 105 model atmospheres with 

T cff = 5400, 5700, 6000, 6300, 6600 K, 
logg = 3.5, 3.9, 4.3, and 

[Fe/H] = -0.8,-1.2,-1.6,-2.0,-2.4,-2.8,-3.2. 
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Fig. B.2. The photometric index of the H/3 line vs. the index 
measured from the profile of H/3 in UVES spectra, l-cr error 
bars: cr(/3(UVES)) = ±0.015 and <j(/3(phot)) = ±0.006 are shown 
in the lower right corner. Stars having [Fe/H] < —2.0 are plotted 
with filled circles, and stars having —2.0 < [Fc/H] < —1.0 with 
open circles. 



Appendix C: S and Zn abundances in DLAs 

Table C.l provides a compilation of S and Zn abundances in 
20 DLA systems as derived from high-resolution spectra ob- 
tained with 8-10 m class telescopes. References to the origi- 
nal works are given in the last column of the Table. S abun- 
dances have been determined from the AA1250, 1253, 1259 
Sn triplet and Zn abundances from the AA2026,2062 Znii 
doublet. Both Sn and Znii are the major ionization stages 
of their respective elements in Hi regions, and corrections 
for unobserved ion stages are expected to be unimportant 
(e.g. Vladilo et al. 2001). Neither S nor Zn show much affin- 
ity for dust and the problem is further lessened in DLAs 
which generally show only mild depletions of even the re- 
fractory elements (Pettini et al. 1997). Hence, the values 
given in Table C.l should reflect the true interstellar abun- 
dances of these two elements in the high redshift galaxies 
giving rise to the damped Lyman-alpha systems. 
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Fig. B.3. The HP2 index of the US line (from Ryan et al. 1999) 
vs. /3(UVES). 



the HP2 index of the BS line (Beers et al. 1999) corre- 
lates extremely well with /3(UVES) for 19 stars in common 
with Ryan et al. (1999). Hence, for the group of metal-poor 
turnoff stars, one may use the HP2 index, which is based 
on medium-resolution spectra, to determine differential val- 
ues of T e ff with the same high precision as obtained from 
/?(UVES). 

The absolute visual magnitude of a star, My (phot), is 
determined from the Stromgren indices {b — y), mi, c\ using 
a calibration derived by Schuster et al. (2006) for a sample 
of stars with accurate (ct(it)/it < 0.15) Hipparcos paral- 
laxes. Reddening corrections were applied if E(b — y) > 0. 
The derived values of My (phot) are given in column 11 of 
Table B.l. In columns 12 and 13, the absolute magnitude 
derived directly from the Hipparcos parallax and its error 
are given for stars with a parallax error less than 30%. 
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Table B.l. V and K a magnitudes, Stromgren indices, and photometric and spectroscopic indices of the strength of the H/3 line. 
The derived interstellar reddening of (b — y) is given in column 10, and the absolute visual magnitude derived from the Stromgren 
photometry is given in col. 11. Columns 12 and 13 give the absolute visual magnitude and its 1-a error derived from the Hipparcos 
parallax if available with an error less than 30%. 



ID 


V 


K s 


b-y 


mi 


Cl 


/3(phot) 


Ref. a) 


/3(UVES) 


E(b - y) 


My (phot) 


My (par) 


a 


BD -13°3442 


10.288 


9.018 


0.308 


0.050 


0.385 


2.622 


4 


8.820 


0.022 


3.67 






CD -30° 18140 


9.949 


8.655 


0.323 


0.047 


0.344 


2.606 


1 


8.755 


0.019 


3.88 


4.19 


0.46 


CD -35° 14849 


10.568 


9.293 


0.321 


0.040 


0.293 


2.603 


1 


8.751 


0.011 


4.47 






CD -42° 14278 


10.216 


8.770 


0.361 


0.040 


0.229 


2.576 


1 


8.587 


0.015 


4.94 






G 11-44 


11.091 


9.740 


0.335 


0.061 


0.248 


2.588 


1 


8.660 


-0.005 


4.80 






G 13-09 


9.998 


8.739 


0.311 


0.048 


0.373 


2.609 


1 


8.806 


0.021 


3.75 


3.70 


0.59 


G 18-39 


10.392 


9.017 


0.346 


0.073 


0.286 


2.581 


1 


8.638 


0.002 


4.35 






G 20-08 


9.948 


8.498 


0.356 


0.047 


0.251 


2.574 


1 


8.671 


0.020 


4.79 


4.47 


0.43 


G 24-03 


10.467 


9.020 


0.363 


0.057 


0.271 


2.585 


1 


8.616 


0.019 


4.53 






G 29-23 


10.230 


8.831 


0.339 


0.059 


0.332 


2.590 


1 


8.706 


0.017 


3.87 






G 53-41 


11.022 


9.595 


0.356 


0.083 


0.271 


2.589 


1 


8.580 


-0.003 


4.54 






G 64-12 


11.459 


10.208 


0.307 


0.043 


0.337 


2.617 


4 


8.860 


0.020 


4.37 






G 64-37 


11.144 


9.923 


0.300 


0.054 


0.333 


2.623 


4 


8.861 


0.007 


4.31 






G 66-30 


11.028 


9.789 


0.305 


0.072 


0.358 


2.634 


1 


8.962 


0.016 


4.20 






G 126-62 


9.478 


8.075 


0.330 


0.063 


0.327 


2.588 


1 


8.739 


0.012 


4.00 


4.06 


0.37 


G 186-26 


10.829 


9.590 


0.306 


0.041 


0.339 


2.608 


1 


8.844 


0.019 


4.31 


5.12 


0.50 


HD 106038 


10.179 


8.761 


0.342 


0.092 


0.264 


2.583 


1 


8.595 


-0.017 


4.69 


4.99 


0.36 


HD 108177 


9.671 


8.354 


0.330 


0.059 


0.287 


2.594 


1 


8.665 


-0.002 


4.36 


4.87 


0.26 


HD 110621 


9.906 


8.566 


0.337 


0.067 


0.324 


2.595 


1 


8.679 


0.007 


3.89 


4.12 


0.46 


HD 140283 


7.213 


5.588 


0.380 


0.033 


0.284 


2.564 


1 


8.434 


0.013 


3.83 


3.36 


0.12 


HD 160617 


8.733 


7.311 


0.347 


0.051 


0.331 


2.584 


1 


8.594 


0.010 


3.51 


3.37 


0.31 


HD 179626 


9.210 


7.680 


0.373 


0.095 


0.293 


2.588 


1 


8.551 


0.005 


4.29 


3.57 


0.39 


HD 181743 


9.687 


8.274 


0.351 


0.052 


0.224 


2.582 


1 


8.567 


-0.001 


4.97 


4.95 


0.34 


HD 188031 


10.148 


8.857 


0.328 


0.058 


0.310 


2.592 


1 


8.730 


0.009 


4.19 






HD 193901 


8.660 


7.144 


0.383 


0.099 


0.221 


2.568 


1 


8.443 


-0.006 


5.16 


5.46 


0.12 


HD 194598 


8.354 


6.982 


0.344 


0.091 


0.269 


2.588 


1 


8.624 


-0.011 


4.63 


4.62 


0.15 


HD 215801 


10.044 


8.642 


0.334 


0.049 


0.330 


2.572 


1 


8.589 


-0.001 


3.47 






LP 815-43 


10.912 


9.650 


0.304 


0.048 


0.382 


2.623 


1 


8.913 


0.031 


4.18 






CD -24° 17504 


12.124 


10.807 


0.322 


0.043 


0.283 


2.609 


4 


8.770 


0.010 


4.60 






CD -71°1234 


10.44 


9.114 










5 


8.790 


0.019 


3.67 






CS 22943-095 


11.762 


10.469 


0.324 


0.045 


0.335 


2.619 


3 


8.804 


0.025 


4.17 






G 04-37 


11.426 


9.974 


0.363 


0.031 


0.306 


2.614 


4 


8.759 


0.053 


4.43 






G 48-29 


10.467 


9.264 


0.298 


0.057 


0.351 


2.620 


4 


8.911 


0.013 


4.29 






G 59-27 


10.892 


9.560 


0.321 


0.058 


0.301 


2.595 


1 


8.750 


0.004 


4.32 






G 126-52 


10.996 


9.719 


0.322 


0.045 


0.347 


2.608 


2 


8.846 


0.031 


4.20 






G 166-54 


11.005 


9.731 


0.324 


0.043 


0.322 


2.619 


4 


8.842 


0.027 


4.42 






HD 84937 


8.332 


7.062 


0.303 


0.056 


0.354 


2.613 


1 


8.821 


0.009 


3.95 


3.77 


0.19 


HD 338529 


9.370 


8.144 


0.308 


0.045 


0.366 


2.600 


1 


8.832 


0.023 


3.98 


3.48 


0.46 


LP 635-14 


11.362 


9.997 


0.347 


0.026 


0.366 


2.611 


4 


8.820 


0.060 


4.02 






LP 651-4 


12.053 


10.842 


0.321 


0.023 


0.340 


2.615 


4 


8.813 


0.033 


4.27 







a) References of V, b - y, mi, ci and /?(phot): (1) Schuster & Nissen (1988), (2) Schuster et al. (1993), (3) Schuster et al. (1996), 
(4) Schuster et al. (2006), (5) Ryan (1989). For all stars, K s has been taken from the 2MASS catalogue (Skrutski et al. 2006). 
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Table C.l. S and Zn abundance measurements in DLAs 



QSO 


~abs 


log iVYH I) 
(cm' 2 ) 


log ATS II) 

(cm' 2 ) 


log JV(Zn II) 
(cm' 2 ) 


fZn/BT 


[S/Znl b 


Ref. c 


Q0000-2620 


3.3901 


21.41 + 0.08 


14.70 ± 0.03 


12.01 + 0.05 


-2.01 


+0.14 


10 


Q0013-004 


1.97296 


20.83 + 0.05 


15.28 + 0.03 


12.74 ± 0.04 


-0.70 


-0.01 


11 


Q0100+130 


2.30903 


21.37 + 0.08 


15.09 + 0.06 


12.47 + 0.01 


-1.51 


+0.07 


2 


Q0201+365 


2.4628 


20.38 


15.29 + 0.01 


12.76 + 0.05 


-0.23 


-0.02 


14, 12 


Q0407-4410 


2.5505 


21.13 + 0.10 


14.82 ± 0.06 


12.44 ± 0.05 


-1.30 


-0.17 


8 


Q0407-4410 


2.5950 


21.09 + 0.10 


15.19 ± 0.05 


12.68 + 0.02 


-1.02 


-0.04 


8 


B0528-250 


2.8120 


21.11 + 0.04 


15.56 + 0.02 


13.27 + 0.03 


-0.45 


-0.26 


1 


Q0551-366 


1.96221 


20.50 + 0.08 


15.38 + 0.11 


13.02 ± 0.05 


-0.09 


-0.19 


6 


FJ081240.6+320808 


2.6263 


21.35 


15.63 ± 0.08 


13.15 + 0.02 


-0.81 


-0.07 


13 


0841+1256 


2.37452 


20.99 ± 0.08 


14.69 ± 0.05 


12.10 + 0.02 


-1.50 


+0.04 


3 


0841+1256 


2.47621 


20.78 ± 0.08 


14.48 ± 0.10 


11.69 + 0.10 


-1.70 


+0.24 


3 


SDSS1116+4118A 


2.9422 


20.28 ± 0.05 


15.01 + 0.10 


12.40 ± 0.33 


-0.49 


+0.06 


5 


LBQS 1210+1731 


1.89177 


20.63 ± 0.08 


14.96 ± 0.02 


12.40 + 0.05 


-0.84 


+0.01 


3 


Q1331+170 


1.77637 


21.14 + 0.08 


15.08 + 0.11 


12.54 + 0.02 


-1.21 


-0.01 


2 


HE2243-6031 


2.33000 


20.67 + 0.02 


14.88 ± 0.01 


12.22 + 0.03 


-1.06 


+0.11 


7 


B2314-409 


1.8573 


20.90 + 0.10 


15.10 ± 0.05 


12.52 + 0.03 


-0.99 


+0.03 


4 


LBQS 2230+0232 


1.86359 


20.83 ± 0.05 


15.29 + 0.06 


12.72 + 0.05 


-0.72 


+0.02 


3 


Q2231-002 


2.06616 


20.53 + 0.08 


15.10 + 0.15 


12.30 ± 0.05 


-0.84 


+0.25 


2 


Q2343+125 


2.4313 


20.35 + 0.05 


14.71 ± 0.08 


12.45 + 0.06 


-0.51 


-0.29 


9 


Q2343-BX415 


2.5739 


20.98 ± 0.05 


15.38 ± 0.03 


12.90 + 0.06 


-0.69 


-0.07 


15 



a [log (Zn/H) DL A - log (Zn/H) ], where log (Zn/H) = -7.39 (Asplund et al. 2005). 
b [log (S/Zn) DL A - log (S/Zn) ], where log (S/Zn) = 2.55 (Asplund et al. 2005). 



c References — 1: Centurion et al. (2003); 2: Dessauges-Zavadsky et al. (2004): 3: Dessauges-Zavadsky et al. (2006,2007); 4: 
Ellison & Lopez (2001); 5: Ellison et al. (2001); 6: Ledoux et al. (2002); 7: Lopez et al. (2002); 8: Lopez & Ellison (2003); 9: 
Lu et al. (1998); 10: Molaro et al. (2000); 11: Petitjean et al. (2002); 12: Prochaska et al. (2007); 13: Prochaska et al. (2003); 
14: Prochaska & Wolfe (1999); 15: Rix et al. (2007). 



